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ABSTRACT 


A  guide  has  been  developed  that  contains 
the  Research  Computation  Center  standards  for 
documenting,  revising,  classifying,  submitting, 
identifying,  obtaining,  and  evaluating  computer 
prc  rams  at  NRL.  It  is  anticipated  that  the 
use  of  such  standards  will  promote  efficient 
organization  of  an  RCC  Program  Library  and  will 
provide  for  maximum  effective  usage  of  this 
program  library  throughout  NRL. 


PROBLEM  STATUS 

A  final  report  on  one  phase  of  the  problem; 
work  is  continuing  on  other  phases. 
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RESEARCH  COMPUTATION  CENTER  PROGRAM  PUBLICATION  GUIDE 


I.  INTRODUCI'ION 

This  guide  contains  the  Research  Computation  enter  (RCC)  standards 
for  documenting,  revising,  classifying,  submitting,  identifying,  obtainin’ , 
and  evaluating  programs  which  exist  or  are  to  exist  in  the  Control  Data 
Corporation  CO-OP  Users'  Group  (CO-OP)  and/or  RCC  Program  Library.  These 
proposed  standards  are  in  part  patterned  after  those  proposed  by  CO-OP. 

Page  numbers  of  this  Guide  consist  of  the  Chapter  Number  followed  uy 
a  sequential  number;  e.g.,  "2-1"  means  Chapter  II,  page  1. 

Chapter  II  contains  information  on  how  to  specify  zero  and  the 
letter  0,  how  to  specify  blanks,  and  when  to  use  upper  case  and  lower  case 
letters  in  documentation. 

Chapter  III  specifies  the  format  for  program  documentation.  Write-ups 
of  major  programs  of  a  permanent  nature  will  be  published  as  NRL  Reports 
and/or  NRL  Computing  References.  Write-ups  of  other  programs  that  are  of 
a  more  temporary  or  changeable  nature  will  be  published  as  NRL  Memorandum 
Reports  and/or  NRL  Computing  Bulletins.  The  decision  as  to  which  category 
a  write-up  belongs  usually  will  be  made  by  the  program  originator  and  the 
RCC.  Programs  will  also  be  submitted  to  CO-OP  upon  request  by  the 
contributor. 


Chapter  IV  describes  the  program  identification  numbering  system  that 
is  used  to  uniquely  identify  each  program  or  subroutine  deck.  Hi  is  nunuer 
will  also  be  used  foi  filing  purpuses  in  the  RCC  Program  library. 

Chapter  V  specifies  standards  for  internal  documentation  -f  programs 
using  comment  or  remark  cards. 

Chapter  VI  defines  the  minimum  equipment  configuration  that  is 
assumed  in  documentation  of  programs. 

Chapter  Vl'I  specifies  the  procedures  for  supplementing  or  revising 
a  program  write-up. 

Chapter  VIII  specifies  the  procedure  for  submitting  a  program  to 
CO-OP  and/or  RCC  for  inclusion  in  the  RCC  Program  Library. 

Chapter  IX  specifies  the  contents  of  a  program  brief  which  is 
necessary  in  u.  ilr  <..u  on  tain  an  RCC  Identification  Number  for  a  program 
deck.  Addendum  h  contains  a  sample  program  brief. 


Chapter  X  specifies  the  procedure  for  obtaining  a  program  which 
exists  in  the  CO-OP  and/or  RCC  Program  Horary. 

•’hap ter  XT  gives  the  suomi 1  tal  form  required  by  CO-OP  when  submitting 
a  program,  addendum  0  contains  a  sample  submittal  form. 

Chapter  XII  gi.es  tne  evaluation  form  required  by  the  RCC  when  using 
a  CQ-OP  and/or  RCC  Library  program. 

Chapt  er  XIII  1  is*  s  and  defines  1  he  classification  codes  used  by 
CO-OP  and  the  RCC. 

Addendum  A  .■  ives  the  write-up  of  a  subroutine  using  the  recommended 
t'o rma t  g i  •- e n  In  Chapter  I'll. 

Addendum  B  contains  a  sample  program  brief  as  described  in  Chapter  IX. 

Addendum  C  cor. 'ains  a  sample  submittal  form  as  described  in  Chapter  XI. 

it  should  be  no* ed  that  throughout  this  ,-uide,  any  reference  tr  the 
RCC  I’r..,  ram  I.i  t-rary  is  no*  a  reference  to  the  computer  system  library. 

(i.e.,  the  on-line  library  file).  ihe  RCC  Program  Library  is  simply  a 
•o'.  !  ed  ion  •>!'  pro grams  maintained  on  file  by  the  RCC,  where  each  pro.  ram 
to  -k  is  uniquely  identified  by  an  RCC  Identifi  at  ion  Number  as  defined 
in  C.,apier  It. 
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III.  FORMAT  FOR  DOCUMENTING  PROGRAMS  AND  SUBROUTINES 


Shown  below  is  the  basic  outline  of  the  format,  to  be  used  whenever 
possible,  to  document  programs  or  subroutines.  Following  the  basic 
outline  is  the  format  including  a  detailed  description  of  each  heading. 
This  format  is  patterned  after  that  contained  in  the  "CO-OP  Handbook, 

Part  III:  Program  Catalog,  Section  1:  Program  Classification  Requirements" 
of  June  1968.  Users  wishing  to  submit  their  programs  to  the  RCC  for 
distribution  to  CO-OP  and/or  to  other  users  are  encouraged  to  follow 
these  conventions. 

All  items  of  the  format  should  be  completed  by  the  author  where 
possible.  All  information  essential  to  the  usage  of  the  routine  must 
be  included.  When  an  item  is  not  applicable,  or  has  not  been  determined, 
this  should  be  explicitly  stated.  If  necessary,  the  user  may  consult 
with  RCC  staff  members  to  determine  how  to  handle  any  items  in  doubt. 
Throughout  the  write-up,  examples  should  be  used  to  clarify  and 
illustrate  usage. 

Addendum  A  gives  the  write-up  of  a  subroutine  using  the  format 
specified  below. 

The  basic  outline  is  as  follows: 

1.0  IDENTIFICATION 

1.1  Title 

1.2  Identification  Name 

1.3  Classification  Code 

1.4  RCC  Identification  Number 

1.5  Entry  Points 

1.6  Programming  Language 

1.7  Computer  and  Configuration 

1.8  Contributor  or  Programmer 

1.9  Contributing  Organization 

1.10  Program  Availability 

1.10.1  Submittal 

1.10.2  On  File 

1.11  Verification 

1.12  Date 

2.0  PURPOSE 

2.1  Description  of  the  Routine 

2.2  Problem  Background 

3.0  USAGE 

3.1  Calling  Sequence  or  Operational  Procedure 

3.2  Arguments,  Parameters,  and/or  Initial  Conditions 

3.3  Space  Required  (Dec Iraki,  and  Octal) 

3.3.1  Unique  Space 

3.3.2  Common  Storage 

3.3.3  Temporary  Storage 

3„4  Messagi  s  and  Instructions  to  the  Operator 
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3*5  Error  Returns,  Messages,  and  Codes 
3*6  Informative  Messages  to  the  User 

3.7  Input 

3.8  Output 

3.9  Formats 

3.10  External  Routines  and  Symbols 

3.11  Timing 

3.12  Accuracy 

3.13  Cautions  to  Users 

3.14  Program  Deck  Structure 

3.15  References  -  Literature  -  Appendices 

4.0  METHOD  OR  ALGORITHM 

5.0  SLOW  CHART  AND/OR  SOURCE  LANGUAGE  LISTING 

6.0  COMPARISON 

7.0  TEST  METHOD  AND  RESULTS 

8.0  REMARKS 


The  format  including  a  description  of  each  item  is  as  follows: 

1.0  IDENTIFICATION 

1.1  Title 

Give  a  descriptive  title  of  the  routine,  limited 
to  approximately  fifty  characters  which  are  included 
in  the  Fortran  character  set,  to  facilitate  indexing 
and  listing. 

1.2  Identification  Name 


This  descriptor  consists  of  the  following  three 
items  of  information,  separated  by  hyphens: 

(a)  the  two  character  classification  code 
as  defined  in  section  1.3* 

(b)  the  CO-OP  code  for  the  contributing 
organization  (NRL). 

(c)  an  alphanumeric  descriptor  of  the  routine 
not  to  exceed  eight  characters.  This 
name  can  not  contain  embedded  blanks  or 
special  characters. 

The  identification  name  is  used  by  CO-OP  and  the 
RCC  for  further  indexing  and  identification. 

For  example,  given  an  NRL  routine  with  a  classification 
code  of  C3  and  a  descriptor  of  BESSEL,  the  identifi¬ 
cation  name  would  be  C3-NRL-3ESSEL. 
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1.3  Classification  Code 

Hive  both  the  code  designation^ )  and  the  verbal 
translation.  The  code  designation  consists  of  a 
letter,  indicating  the  primary  class,  followed  by 
a  digit  indicating  the  sub-class.  The  available 
classification  codes  and  their  descriptions  are 
given  in  Chapter  XIII. 

1 .4  RCC  Identification  Number 

This  eight  digit  code  number  uniquely  identifies 
each  program  or  subroutine  deck.  The  content  and 
construction  of  this  number  is  described  in 
Chapter  IV .  List  each  program  and  subroutine  along 
with  the  identification  number  for  each  when  appli¬ 
cable  . 

1.5  Entry  Points 

(iiv  '  the  alphanumeric  code  name  of  each  entry 
point  in  alphabetical  order.  If  the  name  of  the 
routine  in  which  the  entry  point  occurs  Is  differ¬ 
ent  fran  the  entry  point  name,  then  the  name  of 
the  routine  should  appear  in  parentheses  to  the 
right  of  the  entry  point  name. 

1 .6  Programming  Language 

Indicate  the  language( 3 )  in  which  the  routine  is 
written,  and  specify  tie  type  of  routine,  i.e., 
program,  subroutine,  ot  system.  Also  specify  the 
operating  system  under  which  the  routine  was  written 
and  tested. 

1.7  Computer  and  Configuration 

Specify  the  computer  for  which  the  routine  was 
written  (e,g.  -  CDC  3800).  If  a  special  configura¬ 
tion  or  equipment  complement  different  from  that 
defined  in  Chapter  VI  is  needed,  state  the  particulars. 

1.8  Contributor  or  Programmer 

Specify  the  name,  NRL  code,  branch,  and  division 
of  the  programmer( s )  who  wrote  the  routine.  If 
someone  other  than  the  prcgrammer(s)  given  should 
be  contacted  in  case  of  difficulties  with  the 
routine  include  this  information. 
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1.9 


Contributing  Organization 


Include  the  CC-OP  organization  code  (NRL),  name 
(Naval  Research  Laboratory),  and  address  (Washington, 
D.C.  20390)  of  the  contributing  organisation. 

1.10  Program  Availability 

1.10.1  Submittal 

List  all  of  the  materials  pertaining  to  the 
routine  which  are  being  submitted  to  CO-OP 
and/or  RCC  (Example:  Program  write-up, 

Fortran  source  deck,  source  listing,  etc.). 

1.10.2  On  File 

Specify  means  of  obtaining  the  routine  if 
other  than  through  the  means  stated  in 
Chapter  X.  For  example,  the  routine  might 
be  on  the  on-line  system  library,  in  a  perma¬ 
nent  file  on  mass  storage,  or  on  file  by  the 
contributor.  Also  specify  all  the  materials, 
pertaining  to  the  routine,  which  exist  on  file 
by  the  contributor  if  other  than  those  given 
in  Section  1.10.1. 

1.11  Verification 

this  section  is  to  oe  used  by  the  RCC  to  indicate  to 
what  extent  the  RCC  has  checked  the  routine-  and  the 
testing  procedures  used. 

1.12  Date 

Specify  the  date  of  the  routine  write-up.  Should  any 
revisions  be  made  to  sn  existing  routine,  the  original 
date  and  the  dates  of  any  other  revisions  shall  be 
included  together  with  the  date  of  the  current  revision. 


2.0  PURPOSE 


2.1  Description  of  the  Routine 

Give  a  brief  general  description  of  what  the  routine 
does.  Hie  purpose  should  provide  sufficient  information 
so  that  a  reader  can  decide  whether  or  not  the  routine 
may  meet  a  particular  need. 


. ::  Pn •:  'em  i-aekcrounj 

- St—  — — — 

Specify  the  unique  problem  characrerist ics  that  leu 
!  -•  the  program,  prior  techniques  employed  in  ar.  at  ■ 
to  resolve  the  protlem,  and  1  imitations  of  prior 
techniques  that  prevented  ar,  effective  solution. 

j • c  USAGE 

3 . 1  Calling  Sequence  or  Operational  Procedure 

Specify  the  eall(s;  ”  procedures )  which  is  necessary 
to  utilize  the  reutir.e(s),  including  the  use  of  all 
entry  points.  Specify  and  define  the  function  of  all 
entry  points. 

3.2  Arguments,  parameters,  and/or  Initial  Conditions 

Describe  all  arguments,  parameters,  and  arrays  needed 
to  utilize  the  reutine(s).  This  description  should 
include  what  the  arguments  represent  in  the  physical 
problem,  mode,  and  any  size  information,  or  limitations. 
Also  describe  any  assumptions  made  by  each  routine 
wnen  it  is  called. 

3 • 3  Spa  -b  Required  (Decimal  and  Octal) 

3.3.1  Unique  Storage 

Specify  all  unique  core  space  requirements  of 
the  routine,  excluding  common  blocks. 

3 -3.2  Common  Blocks 

Specify  the  name  and  length  of  all  common  Hecks 
required  by  the  routine.  Also  specify  the  names 
of  the  arrays  or  cells  within  each  common  b’.ck, 
including  their  length,  their  function,  their 
mode,  and  all  routines  which  use  each  cell  or 
irray.  One  concise  way  to  show  some  of  this 
information  is  shown  by  the  following  imaginary 
exampl e : 
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ROUTINES  COMMON  SLOCKS 


— 

BLOCK  RDC 

BLOCK  SRI 

1 

ADYG 

ADV0 

* 

NCW 

NTP 

# 

S0RTAR(257) 

* 

scale( 96 ) 

1 

ADDLIST 

X 

X 

X 

y-1 - 

i  FI 

1 

X 

— 

X 

MUI.C 

X 

X 

X 

X 

j  RDSPEC 

X 

X 

j  SPECPR 

i 

X 

X 

X 

X 

!  SPECPN 

X 

X 

r 

TOTAL 

X 

X 

*  indicates  real  (floating  point)  node 
(all  other  variables  are  in  integer  mode) 


3-3-3  Temporary  Storage 

Temporary  storage  areas  are  those  areas  that 
contain  pertinent  information  only  at  certain 
times  in  a  run  and  at  all  other  times  the  contents 
of  the  areas  could  be  destroyed  without  consequence. 
Specify  the  names  of  temporary  storage  and  the  size 
of  each  along  with  the  names  of  all  routines  that 
use  the  temporary  storage.  It  should  be  approp¬ 
riately  noted  if  the  temporary  storage  is  included 
in  unique  or  common. 

3.4  Messages  and  Instructions  to  the  Operator 

Specify  and  define  any  messages  given  to  the  operator  by 
the  routine  and  any  special  operator  instructions  such  as 
special  tape  handling  procedures.  In  general,  output 
messages  to  the  operator  are  discouraged  by  the  RCC, 

3.5  Error  Returns,  Messages,  and  Codes 

Include  a  list  of  the  error  situations  (error  returns, 
error  codes,  and  error  stops)  that  are  detected  and  the 
action  taken.  Any  diagnostic  messages  should  be  speci¬ 
fied  and  described.  Also  Bpecify  whether  each  error 
causes  normal  or  abnormal  termination  of  job  execution. 
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If  an  error  situation  in  this  routine  can  cause 
diagnostic  action  by  another  routine  or  the  system 
monitor,  this  action  should  be  described.  When 
possible,  recommend  the  action  for  the  programmer 
and/or  operator  to  take. 


Informative  Messages  to  the  User 


Specify  and  define  any  messages  that  give  information 
rather  than  indicate  an  error.  All  messages  in  this 
category  assume  that  the  run  continues  without  any 
action  on  the  part  of  the  user. 


3.7  Inpn  t 

Include  (where  applicable)  the  input  media  (cards, 
magnetic  tape,  etc.),  logical  unit  assignments, 
labels,  densities,  tape  numbers,  and  clarifying 
descriptions  for  all  input  required  by  the  routine. 

3.8  Output 


Include  (where  applicable)  the  output  media,  logical 
unit  assignments,  labels,  densities,  tape  numbers, 
and  clarifying  descriptions  for  all  output  generated 
by  the  routine.  Estimates  of  the  output  volume  should 
also  be  included  when  possible. 


3.9  Formats 

Describe  the  formate  and  contents  of  all  input  and 
output  data.  Often  it  is  helpful  to  include  a  list 
of  format  statements  or  format  specifications  for 
each  field. 


External  Routines  and  S\ 


List  in  alphabetic  order  each  routine  and  the  external 
symbols  that  each  routine  calls.  One  concise  way  to 
show  this  information  is  by  use  of  a  table  that  would 
have  the  names  of  the  routines  down  the  left  side,  all 
external  symbols  across  the  top  and  checks  in  the 
appropriate  rows  and  columns.  If  the  name  of  the 
routine  in  which  the  external  symbol  occurs  is  differ¬ 
ent  from  the  external  symbol,  then  the  name  of  the 
routine  should  appear  in  parentheses  to  the  right  of 
the  external  symbol.  Also  indicate  where  each  external 
routine  is  physically  located  (system  library,  RCC 
Program  Library,  etc.)  if  not  within  the  program  deck. 
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j.U  Timing 

When  possible,  estimates  of  running  time  should  le 
pro . ided. 

3.12  Accuracy 

Per  mathematical  routines,  statements  of  accuracy, 
precision,  or  the  magnitude  and  propagation  of  errors 
should  be  included. 

3.13  Cautions  to  Users 

Include  anything  of  a  cautionary  nature. 

3.1 4  Program  Deck  Stricture 

Specify  the  decjc  structure  for  a  run,  including  control 
cards,  routines,  and  data. 

3-15  References  -  Literature  -  Appendices 

List  and  describe  any  references,  literature,  appendices, 
or  other  documentation  that  pertain  to  this  ro^tiiie. 


4.0  METHOD  OR  ALGORITHM 

Describe  how  the  routine  works  in  regard  to  the  mathematics 
or  algorithm.  Describe  the  programming  teehr.iq  ;es  and  methods, 
supporting  theory,  design,  computation  equations,  and  the 
derivations  of  the  computation  equations  to  substantiate  or 
illustrate  the  program  and  emphasize  its  advantages,  if  any, 
over  previous  methods. 


5.0  FLOW  CHART  AND/OR  SOURCE  LANGUAGE  LISTING 

Flow  charts  should  be  submitted  whenever  they  are  avail sole 
and  not  too  lengthy,  particularly  for  machine  language  routines. 
Large  programs  should  nave  both  general  and  detailed  flow 
charts.  The  general  flow  chart  consists  of  one  block  for 
each  subroutine.  A  source  language  listing  should  only  be 
included  when  it  is  fairly  short.  If  a  flow  chart  or  source 
language  listing  is  available  but  is  not  being  included  here, 
state  how  it  may  be  obtained. 
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_  .  RIP  OH 

if  similar  programs  exist  in  the  CO-OP  and/or  RCC  program 
1C : rarv  a  comparison  should  be  made  with  respect  to  such 
programs  discussing  the  advantages  and  disadvantages  of  the 
t  it  *w  routine  and  indicating  any  programs  whi  di  are  made 

o :  u.  !  e '  e . 


7 . 0  Ttuii  MKTHOP  AND  RESULTS 

In  -I ude  test  data  used  and  results  obtained  in  checking  tne 
.  peratlon  and  accuracy  of  tne  routine. 

b .  c  hKKRRKO 

dive  any  older  information  that  may  be  useful. 


3-9 


IV.  PROGRAM  IDENTIFICATION  NUMBERING  SYSTEM 


An  RCC  program  identification  number  is  given  to  any  program  or 
subroutine  neck  *.;•  be  filed  in  the  RCC  Program  Library.  This  number 
provides  a  means  for  identifying  the  following  items: 

1 .  Source  program  or  subroutine  decks 

2.  Application  classification  of  the  program 

3.  The  program  originator 

4.  The  level  of  program  revision 

Ihe  RCC  program  identification  number  consists  of  eight  characters 
that  are  determined  as  follows: 


Character 

Number 

Function 

Current 

Assignments 

1,2 

Classification 

Code 

One  letter  followed 
by  one  digit  (Sec  tii 
in  Chapter  XIII ) 

A  3-digit  serial 
(within  each 
class)  number 
assigned  oy  RCC 

000-999  (range) 

6 

Originator  code 

C  -  CO-OP 

R  -  RCC 

0  -  Other  NRL 

T  -  Tenant  Acti.ities 

N  -  Non -NRL 

7,S 

A  2-digit,  program 
revision  code 

00-99  (range) 

A  program  identification  number  can  be  obtained  from  the  RCC  upon 
submittal  of  a  Program  brief  (see  Chapter  IX )  to  the  RCC  by  the  programmer 
If  4 he  program  deck  is  to  be  physically  located  in  the  RCC  Pr  gram  Library 
file  it  is  additionally  required  that  the  program  card  deck  which  is 
submit. '  ed  must 

(a)  be  properly  sequence -numbered  in  card  columns  7 j-8c  wi*h 
sequential  numbers  increasing  in  value  through  the  deck,  anu 

(b)  cun  bain  internal  documentation  as  given  in  Chapter  \. 

'ihe  ICC  program  identification  number  should  be  included  cit  a  comment 
or  remark  curd  following  the  routine  name,  at  the  beginning  of  t tie  sytruoli* 
deck  (see  Chapter  V).  ihe  identification  number  should  also  be  included 
in  any  documentation  about  the  routine. 


I  he  program  revision  code  is  used  as  follows:  fhe  initial  vers 
of  a  new  program  or  subroutine  deck  is  assigned  a  revision  code  of 
Each  time  the  program  or  subroutine  is  revised  or  modified  in  any  wa 

number  is  .incremented  i  y  one. 


V.  PROGRAM  INTERNAL  DOCUMENTATION  STANDARDS 


Each  program  or  subroutine  deck  should  contain  internal  documentation 
in  tiie  form  of  comment  or  remark  cards.  This  internal  documentation  should 
perform  three  functions : 

1.  Identify 

2.  Define  usage 

j.  Clarify  program 

Identifying  comments  or  remarks  are  strongly  recommended  for  all 
programs  and  are  required  by  the  RCC  before  a  routine  will  be  put  in  : he 
RCC  Program  Library  or  accepted  for  publication.  These  comments  or  remarks 
should  appear  immediately  following  the  routine  name  and  should  con'air.  *  .he 
following  information  in  the  specified  order: 

IDENT  NUMBER  - 
TITLE  - 
IDENT  NAME  - 
LANGUAGE  - 
COMPUTER  - 
CONTRIBUTOR  - 
ORGANIZATION  - 
DATE  - 
PURPOSE  - 

The  above  headings  are  defined  in  Chapter  III,  Sections  1.0  and  2.0.  File 
information  supplied  in  these  comments  or  remarks  need  not  be  as  detailed 
as  in  the  write-up.  The  above  identifying  headings  should  be  included  as 
part  of  the  comments  or  remarks. 

Usage  information  comments  or  remarks  are  not  required  cut  are  s*r.  ugly 
encouraged,  'these  comments  or  remarks  should  follow  the  identifying  -\mments 
or  remarks  and  should  include  the  following  type  of  information: 

CALLING  SEQUENCE  - 
ARGUMENTS  - 
I/O  INFORMATION  - 
EXTERNAL  SYMBOLS  - 

The  above  headings  are  defined  in  Chapter  III,  Section  3«0*  The  infont.u'  ion 
supplied  in  these  comments  or  remarks  need  not  be  as  detailed  as  in  the 
write-up.  'Ihe  above  headings  should  be  included  as  part  of  the  cornmer. ' s  or 
remarks . 

Clarifying  comments  or  remarks  should  appear  throughout  t lie  rc.  i.ue 
to  describe*  what  the  routine  is  doing  at  any  gi-en  point. 

In  addition  to  the  above  documentation,  it  is  required  that  each,  deck 
of  a  routine  be  properly  sequence-numbered  in  card  columns  73-^C  with  v.e 
numbers  increasing  in  value  through  the  deck. 
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VI .  MINIMUM  EQUIPMENT  CONFIGURATION 


The  minimum  equipment  configuration  assumed  for  all  documentation 
routines  is  the  basic  CDC  38OO  computer  with  32K  core  storage,  2  magnet 
tape  units,  drums,  card  punch,  card  reader,  and  a  line  printer.  If  a 
routine  requires  equipment  in  addition  to  that  stated  above,  such  shoul 
be  stated  in  the  documentation  (Chapter  III,  Section  1.7).  The  above 
minimum  configuration  is  similar  to  that  given  by  CO-OP. 
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VII.  PROGRAM  ADDENDA  AND  REVISIONS 


When  an  error  is  found  an  RCC/CO-OP  Library  Evaluation  Form  (see 
Chapter  XI I )  should  be  sent  to  the  RCC  Program  Library,  Code  781 3L. 

This  evaluation  should  explain  the  error  in  detail.  To  correct  or  add 
to  an  existing  routine,  use  the  documentation  format  given  in  Chapter  III 
and  follow  the  standard  program  submittal  procedure  described  in  Chap¬ 
ter  VIII.  Only  those  items  that  are  affected  need  be  resubmitted.  Only 
the  originator  (unless  he  is  unavailable)  of  a  program  can  make  corrections 
or  modifications  to  the  CO-OP  and/or  RCC  Program  Library  file  copy  of  a 
program. 

An  addendum  to  a  publication  will  be  indicated  by  a  letter  following 
the  original  publication  number  (e.g.  -  A,  B,  etc.,  for  each  succeeding 
addendum). 

A  revision  of  a  publication  will  normally  be  accomplished  by 
re-publishing  the  entire  document.  A  revision  of  a  publication  will  be 
indicated  by  .n  following  the  RCC  publication  number,  where  n  =  1,2,.... 
for  each  succeeding  revision.  However,  a  revision  to  an  NRL  Report  or  an 
NKL  Memorandum  Report  will  be  given  a  new  report  number,  with  the  old 
number  specified  as  being  superseded.  A.  revision  should  specify  the  dates 
of  the  original  and  all  revisions. 


VIII. 


PROGRAM  SUBMITTAL  PROCEDURE 


Users  having  programs  cf  general  interest  which  have  been  thoroughly 
tested  and  documented  may  submit  the  programs  to  the  RCC  for  publication 
as  an  NKL  Computing  Note,  Bulletin,  or  Reference  (Programs  to  be  sent  to 
CO-OP  or  to  be  documented  as  NRL  Reports  or  Memorandum  Reports  must  follow 
additional  standard  procedures.)*  I'16  programs  should  be  written  in  a 
compiler  or  assembly  language  generally  available  through  the  RCC.  A  memo¬ 
randum  should  be  written  to  the  RCC  Program  Library,  Code  T813L  requesting 
that  the  program  be  published  and  included  in  the  RCC  Program  Library. 
Included  with  the  memorandum  should  be  the  following: 

1.  Program  Brief  (defined  in  Chapter  IX;  example  in  Addendum  B) 
if  the  KCC  does  not  already  have  one  for  the  routine(s).  If  a 
Program  Brief  is  currently  in  existence,  the  RCC  Identification 
Number  for  the  routine  must  be  supplied. 

2.  Program  Source  Listing  resulting  from  the  assembly  or  compilation 
of  the  submitted  source  deck. 

3.  Program  Source  Deck  on  cards,  interpreted,  sequence-numbered,  and 
internally  documented  (see  Chapter  V);  or  the  RCC  Identification 
Number  of  the  routine(s)  if  already  in  existence  in  the  RCC  Program 
Library  file. 

4.  Program  write-up  preferably  in  the  form  specified  in  Chapter  III. 
See  the  example  in  Addendum  A. 

If  a  program  is  to  be  sent  to  CO-OP,  the  CO-OP  Submittal/Revision  Form 
must  also  be  supplied  (see  Chapter  XI  and  Addendum  C). 


IX.  PROGRAM  BRIEF 


In  order  to  obtain  an  RCC  Identification  Number  for  a  program  or 
subroutine  d ck ,  a  Program  Brief  must  be  supplied  to  the  RCC.  The  Program 
Brief  shoula  contain  the  following  information: 

TITLE. 

IDENTIFICATION  NAME. 

PROGRAMMING  LANGUAGE. 

COMPUTER  AND  CONFIGURATION. 

PROGRAMMER. 

NHL  ORGANIZATION. 

DATS. 

PURPOSE. 

PROGRAM  STATUS. 

The  above  headings  are  described  in  part  by  Chapter  III,  Sections  1.0 
and  2.0.  The  NRL  ORGANIZATION  section  should  include  the  NRL  Division, 
Branch,  and  Code,  in  that  order.  The  PROGRAM  STATUS  section  should 
indicate  what  the  current  status  of  the  program  is  (complete,  incomplete, 
etc. ). 

A  sample  Program  Brief  has  been  prepared  in  Addendum  B. 


X.  PROCEDURE  FOR  OBTAINING  A  PROGRAM 

All  programs,  for  which  documentation  exists  in  the  form  specified 
in  Chapter  III,  will  be  available  in  the  RCC  Program  Library  unless  stated 
otherwise  in  the  write-up  (Chapter  III,  Section  1.10,2).  To  obtain  a 
program  existing  in  the  CO-OP  and/or  RCC  Program  Library  contact  Research 
Computation  Center  Program  Library  personnel.  Users  of  all  routines 
obtained  from  the  CO-OP  and/or  RCC  Program  Library  are  requested  to  submit 
an  RCC/C0-0P  Library  Evaluation  Form  (see  Chapter  XIII )  to  the  RCC  Program 
Library,  Code  7813L.  This  evaluation  should  contain  detailed  usage  informa¬ 
tion  about  discovered  errors  or  inconsistencies,  and  any  other  pertinent 
information. 
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XI .  USER  GROUP  LIBRARY  SUBMITTAL /REVISION  FORM 

The  Piv-'ram  Abstract  Cover  Sheet  as  given  at  the  end  of  this  chapter 
must  be  filled  out  whenever  a  program,  program  revision,  or  write-up 
revision  is  neing  submitted  tc  CO-OP.  A  sample  Program  Abstract  Cover  Sheet 
has  been  prepared  in  Addendum  C.  The  categories  on  the  form  are  defined 
as  follows: 


1.  Library:  Th*  ^ord  "CO-OP"  should  be  placed  in  the  blank 
preceding  the  word  "LIBRARY"  to  indicate  that  the  routine  should  be  place! 
in  the  CO-OP  Library  as  opposed  to  SWAP,  POOL,  etc. 

2.  Catalog  Identification:  The  "Identification  Name"  information 
as  described  in  Chapter  III,  Section  1.2. 

3.  Title  (Descriptive):  See  Chapter  III,  Section  1.1. 

4.  Source  Language(s)  or  Dialect(s);  Name  of  specific  source 
language(s)  and  dialect(s),  or  versions,  in  which  the  program  is  coded. 

5.  Operating  System  and  Version:  Name  of  Operating  System  (e.g., 
DRUM  S COPE;  and  Version  (e.g.,  Version  2.0). 

6.  Computer  and  Configuration:  The  computer  and  configuration  for 
which  the  program  has  been  written  and  checked  out  should  be  given.  This 
should  include  the  storage  space  required. 

7.  Contributing  Organization:  The  full  name  of  the  contributing 
organization  as  defined  in  Chapter  III,  Section  1,9. 

8.  Submitter/Programmer:  The  name(s)  of  the  programmer ( s )  who 
wrote  the  program.  Include  the  programmers )  NHL  Code(s). 

9*  Revisor:  The  name(s)  of  the  programmers )  who  revised  the 
program.  This  is  left  blank  when  the  original  program  is  being  submitted. 

10.  Date  Written:  The  date  of  the  original  program  shall  be 
specified.  When  submitting  a  revision,  the  date  of  the  revision  should 
also  be  given. 

11.  Program  Material :  The  number  of  pages  in  the  program  write-up 
(including  the  Program  Abstract  Cover  Sheet)  should  follow  "Write-up 
Page  Count". 

If  a  card  deck  is  submitted,  the  number  of  cards  in  the  deck  should 
follow  "Sequenced  Source-Card  Count".  If  the  program  is  submitted  on 
magnetic  tape  it  should  consist  of  one  physical,  file  and  the  total  number 
of  BCD  card  images  should  be  entered  following  "M/T  Record  Count". 

If  the  program  is  submitted  on  paper  tape,  the  number  of  recognizable 
records  on  the  tape  should  be  entered  following  "P/T  Record  Count". 
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12.  Qt-iier  Available  Documents :  The  submitter  should  specify 
when  other  items  are  submitted,  such  as  listings,  flow-charts,  sample 
print-outs,  etc.  A  listing  need  not  be  included  if  the  card  deck  is 
sequenced. 

13.  Original /Revised  Program  Catalog  Abstract:  The  submitter 
should  type  his  aisiraet  exactly  as  he  wishes  it  to  appear  in  forthcoming 
issues  of  user  group  program  catalogs.  The  abstract  should  be  no  lon.-er 
than  two  short  paragraphs  and  should  contain  the  following  information: 

Purpose 

Entry  Points  (in  alphabetical  order) 

Program  Type  (subroutine,  complete  routine,  o’*  system) 

14.  Nature  of  tie  vision  (or  Additional  Information):  This  space  *s 
provided  for  giving  more  detailed  information  about  a  revision  than  is 
needed  for  the  abstract  in  the  program  catalog. 

All  information  provided  on  the  Program  Abstract  Cover  Sheet  must 
be  in  standard  48-character  code,  for  keypunching.  Therefore,  mathematical 
symbols  and  other  special  characters,  though  permissible  elsewhere  in  the 
aocumentation,  must  be  converted  ’or  typing  on  the  Program  Abstract  Cover 
Sheet. 

Copies  of  the  Program  Abstract  Cover  Sheet  are  available  from  the  RCC. 
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LIBRARY 


PROGRAM  ABSTRACT  COVER  SHEET 


(SUBMITTAL  P^RM) 


CATALOG  IDENTIFICATION 


REVISION 


TITLE  (DESCRIPTIVE) 


SOURCE  LANGUAGE (S) 
AND  DIALE CT(S)  : 

OPERATING  SYSTEM 
AND  VERSION  : 


COMPUTER  A  CONFIGURATION 


CONTRIBUTING  ORGANIZATION  : 


SUBMITTOR/PROG RAMMER  : 


REVISOR  : 


DATE  WRITTEN  : 


ORIGINAL 


REVISION 


PROGRAM  MATERIAL  :  WRITE-UP  PACE  COUNT  s 


SEQUENCED  SOURCE -CARD  COUNT  : 
M/T  RECORD  COUNT  :  _ 


P/T  RECORD  COUNT  s 


OTHER  AVAILABLE  DOCUMENTATION 


For  Coaplttion  by  Dlatrlbutlon  Aiancy  : 
Original /Pravlou*  Ravi* Ion  <  UGN  _ 
Latast  Ravlalon  :  UGN 
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XII.  RCC/CO-OP  LIBRARY  EVALUATION  FORM 

Each  user  of  an  RCC  or  CO-OP  Library  program  is  requested  to  submit 
an  RCC/CO-OP  Library  Evaluation  Form,  as  given  at  the  end  of  this  chapter, 
to  the  RCC  Program  Library,  Code  7&L3L*  This  evaluation  should  specify 
in  do ‘ail,  any  difficulties  encountered  and  any  corrective  measures  taken 
«h.~e  using  the  program. 
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RCC/CO-OP  LIBRARY  EVALUATION 


DATE 

RCC/CO-OP  IDBJTI  FI  CATION _ 

PROGRAM  TITLE _ 

LANGUAGE _ _ 

usnt's  application* 


DIFFICULTIES  ENCOUNTERED,  IF  ANY* 


COMMENTS,  IF  ANT* 


USER'S  NAME _ 

CODE _ 

TEL.  EXT _ 

PLEASE  SBJD  COMPLETED  FORM  TO  RCC  PROGRAM  LIBRARY,  CODE  70131. 
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XIII.  CLASSIFI CATION  CODES 


The  following  classification  system  is  used  by  the  RCC  to  categorize 
computer  subroutines  and  programs.  This  classification  system  is  almost 
identical  to  that  currently  proposed  by  CO-OP  in  the  "Users  Organizations 
for  Control  Data  Computer  Systems  Handbook,  Part  III:  Program  Catalog, 
Section  1:  Program  Classification  Requirements",  June  1968.  The  classifi¬ 
cation  code  consists  of  a  letter,  indicating  the  primary  class,  followed 
by  a  digit  indicating  the  sub-class.  Under  each  primary  class  will  be  a 
zero  code  having  the  heading  "General"  under  which  will  be  filed  the 
routines  that  are  not  appropriate  for  the  listed  secondary  classes. 

The  primary  classes  are  listed  below  after  which  the  entire 
classification  system  is  given  including  a  description  of  some  of  the 
primary  and  sub-classes. 

A.  ARITHMETIC  ROUTINES 

B.  ELEMENTARY  FUNCTIONS 

C.  POLYNOMIALS  AND  SPECIAL  FUNCTIONS 

D.  OPERATIONS  ON  FUNCTIONS  AND  SOLUTIONS  OF  DIFFERENTIAL  EQUATIONS 

E.  INTERPOLATION  AND  APPROXIMATIONS 

F.  OPERATIONS  ON  MATRICES,  VECTORS  .AND  SIMULTANEOUS  LINEAR  EQUATIONS 

G.  STATISTICAL  ANALYSIS  AND  PROBABILITY 

H.  OPERATIONS  RESEARCH  TECHNIQUES,  SIMULATION  AND  MANAGEMENT  SCIENCE 

I .  INPUT 

J.  OUTPUT 

K.  INFORMATION  TRANSFER 

L.  EXECUTIVE  ROUTINES 

M.  DATA  HANDLING 

N.  DEBUGGING 

O.  SIMULATION/EMULATION  OF  COMPUmERS  AND  DATA  PROCESSING  INTERPRETERS 

P.  DIAGNOSTICS 

Q.  S'CRVI  ;e  OR  HOUSEKEEPING:  PROGRAMMING  AIDS 

R.  NON-NUMERICAL  MATHEMATICS 

S.  INFORMATION  RETRIEVAL 

T.  APPLICATIONS  AND  APPLICATION-ORIENTED  PROGRAMS 

’J.  LINGUISTICS  AND  LANGUAGES 

V.  GENERAL  PURPOSE  UTILITY  SUBROUTINES 

X.  DATA  REDUCTION 

Y.  USER  ASSIGNED 

Z.  ALL  OTHERS 


The  entire  classification  system  is  as  follows: 


A.  ARITHMETIC  ROUTINES  Routines  of  the  form  of  sub- 

routines  or  partial  programs 
to  be  incorporated  into  programs 
or  subroutines  for  the  purpose 
of  doing  arithmetic  operations 
within  a  given  algebraic  stricture. 

AO.  General 


A1 .  Real 


A2.  Complex 


A3 .  Byte 


Ah.  I/O  Routines 


B.  ELEMENTARY  FUNCTIONS 


Routines  which  perform  exclusively 
real  arithmetic  operations, 
including  single  precision, 
extended  precision,  and  multiple 
precision  representations  of 
real  numbers  (real  opposed  to 
imaginary). 

Routines  which  perform  complex 
arithmetic  operations  including 
extended  and  multiple  precision 
representations  on  complex 
numbers . 

Routines  which  perform  single 
or  multiple  precision  byte 
arithmetic. 

I/O  Routines  designed  for  use 
with  (e.g.)  multiple  precision 
arithmetic  packates  frequently 
contain  entries  for  noth  input 
and  output,  hence  could  legiti¬ 
mately  be  filed  under  I  or  J. 
However,  their  real  intent  and 
purpose  is  for  use  only  with 
special  routines  already  files 
under  A  (filing  under  I*  or 
therefore  merely  loses  them  for 
prospective  users).  While  -his 
type  of  i/o  package  is  not  pro¬ 
grammed  arithmetic  in  -he  -rue 
sense  of  the  term,  its  intimate 
relationship  with  A-eatevor y 
routines  merits  that  classifi  -at  ion. 

Elementary  functions  are  defined 
to  be  those  generated  cy  constants 
and  the  independent  variable  iy 
means  of  a  finite  sequence  of 
elementary  operations.  El  omen: ary 
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operations  on  functions  F(X)  and 
G(X)  are  those  that  yield  any  of 
the  following  F(X)  =  G(X),  F(X)* 
G(X),  F(X)/G(X),  (F(X))**A, 

A**(x),  LOG  (X),  T  (X)  where  T 
is  any  trigonometric  or  inverse 
trigonometric  function  and  X  is 
any  number  (fixed-point,  floating¬ 
point,  complex,  multiple-precision, 
etc. ). 

BO.  General 

Bl.  Trigonometric  and 

Inverse  Trigonometric 

B2.  Hyperbolic  and  Inverse 
Hyperbolic 

B3.  Exponential  and  Logarithmetic 

Refers  to  roots  of  quantities 
not  polynomials. 

Routines  which  include  operations 
on  and  evaluation  of  polynomials 
ar.d  functions  often  referred  to 
as  special  functions,  such  as 
hypergeometric,  gamma,  Bessel , 
elliptic,  etc.  Simultaneous 
nonlinear  equations  shall  include 
operations  on  and  solutions  of 
systems  of  nonlinear  algebraic  and 
transcendental  equations  but  shall 
exclude  differential  type  equations. 

CO.  General 

Cl.  Evaluation  of  Polynomials 

CP.  Roots  of  Polynomials 

C3.  Evaluation  of  Special  Functions 

CL.  Simul taneoua  Nonlinear  Algebraic  Equations 

C<y.  Simultaneous  Transcendental  Equations 


B4.  Roots  and  Powers 


C.  POLYNOMIALS  AND  SPECIAL 
FUNCTIONS 


lo-i 


OPERA I I QNS  ON  FUNCTIONS  AND  SOLUTIONS  OF  DIFFERENT IAL  BQJA'ICKc 

This  a* e gory  shall  incl..do 
routines  to  perf  rm  nutr.eri.-al 
integration  and  r.umeri -al  differ¬ 
entiation  as  well  as  routines  to 
solve  differentia;  equations  loth 
ordinary  and  partial.  Other 
routines  whi ?h  e  alua’e  .  rera*  ions 
of  functions  or  per  fore.  *  he  o rer¬ 
an  ot.  on  a  f motion  shall  air-  tr¬ 
ine!  tided . 


DO.  General 

Dl .  Numerical  Integration 

D 2.  Numerical  Solutions  to 
Ordinary  Differential 
Equations 

Dj.  Numerical  Solutions  to  Partial 
Pi fferential  Equations 

Pi*.  Numerical  Differentiation 

Interpolation  shall  ir.'l  Me 
routines  which,  whet,  yi  .-r.  a  se* 
of  points,  approximate  alues 
based  or.  polynomial s ,  ratios  of 
polynomial s,  sines  atvdA-r  -v sines, 
exponential  functions,  et  •. 
Approximations  shall  me  a-  routines 
to  pro*',  ide  functions  or  -ur.es  *> 
fit  a  set  of  points  »  r  da*a  : y 
polynomial s,  exponentials,  sines, 
cosines,  etc. 

1-X  .  General 

Ki  .  "a;  if-  Look-up  and 

In**- rpol at  ion 

EG.  .r  c  Fitting  (See  also  Go) 

K 3 .  Go h> tiling  (see  al so  Z2 ) 


I  i-1* 


EU.  Extreme  Values j  Includes  variational  problems. 

Minimizing  or  Maximizing 
a  Fiinction 


F.  OPERATIONS  ON  MATRICES,  VEC- 
TORS  ANB  SIMULTANEOUS  LINEAR 
EQUATIONS 


This  category  shall  include  all 
routines  whose  primary  purpose 
is  the  execution  of  matrix  opera¬ 
tions,  such  as  addition,  multipli¬ 
cation,  inversion  and  transposition. 
Evaluation  of  the  determinant  and 
finding  the  eigenvectors  and  eigen¬ 
values  of  a  matrix  shall  be 
included,  as  well  as  solutions  of 
a  system  of  simultaneous  linear 
equations. 


FO.  ueneral 

Fl .  Matrix  Operations 

F2.  Eigenvalues  and  Eigenvectors 


F'3 .  Determinants 

eh.  Simultaneous  Linear  Equations 


(I.  STATISTICAL  ANALYSIS  AND 
FROBAbILI  n 


This  category  is  meant  to  v  er 
routines  which  are  rasicaliy 
statistical  in  nature,  such  as 
finding  means,  standard  delations, 
etc.,  and  routines  like  curve 
fitting  where  the  routine  is 
programmed  explicitly  for  statis¬ 
tical  purposes. 


10.  ><neral 

>1  .  Ta*a  Re  lu  •*  ion 


Correia*  ion  *2.  ’•  R«>rcsu- 
;  u.  Ai.a.  yais» 

{  See  at  3,  v 

Seq:j.»ft •  ‘  a!  Anal .. v- 1  a 

Analysis  of  ’arlar.  e 
ir,’*1  Series 


Refers  to  the  calculation  of  the 
more  common  statistical  parameters 
su -h  as  mean,  median,  s tan  lari 
!'• .  la'  ion,  et*-. 

Include-#  cur  .e  fi  t  tun,  Vi: i  i:  i:; 

**ni  1  -i  ?  iy  for  a*n?  iat  icai  purpose: 


M  :  te  ‘art  (S«**»  at s.  H)  Ir-clu-Vs  Fandom  N.sr.trr  Enera'ord. 
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H.  OPERATIONS  RESEARCH  TECHNIQUES, 
SIMULATION  AND  MANAGEMENT' 
SCIENCE 


Ho,  General 

Hi  .  Linear  Profrananir.j 

HP.  Non 'Linear  Pr\  orannrvino 

Hi.  Transportation  and 
Network  Codes 

Hk.  Simulation  H-'uel i:*- 


Opera*  ions  resear  *h  will  m-.a*  •  r.  ■ 
application  o‘‘  s-'Ier.ti  fi  •  mv'hocis, 
te  -hniques  ar.i  *ools  *o  *  lit*  sol  ;*  ion 
of  pro:  ler.3  invol  in-*  operations 
of  a  system  so  as  *o  provi  ie  t hose 
in  -on  rol  of  the  system  wi  *h 
opt  itnur.  sol u*  ions  to  the  iro  lem. 
Linear  proe"rajr,n:ir.  1.  a  ot  of 
sol'  ing  proVlems  vi*t  "i  es 

and  constraints  stated  a  linear 
manner.  Simulation  is  * r.e  -hr: quo 
of  desorillnr  the  operational 
systems  ty  means  of  mathoma* i  -al 
models.  Management  sci--r.  :e,  ar.u 
interrelated  areas,  will  re  included 
and  will  cover  such  areas  as 
solutions  to  problems  of  resource 
al location.  It  is  visualized  that 
this  category  (h)  will  cover  such 
things  as  production  and  inventory 
scheduling,  logistics,  queuino 
theory,  renewal  theory,  dynamic 
programming,  itame  theory  and  he  er¬ 
istic  projcramiinn.  Mu.nte  Carl, 
methods  used  in  this  type  of  problem 
will  he  included  al  though  Monte  Carl, 
methods  may  also  appear  in  the 
statistical  cate.vry  (  ''  depend  in.-.: 
upon  *he:r  use. 


P  i  nd  i : .  o  the  test  sol  at  ior.  rc  r. 
amon*  all  solutions  of  a  v;,  stem 
of  1  i near  i nequa!  1  -  i ••  s . 

Sol  •  inp  const rai  tied  p  •  Is.  Isa*  ivr; 
pro:  U-ms  ex  *«rp*  *  n  -s-  wiser**  ’  ne 
obje  •;  i  .e  funcM-n  an: 
constraints  are  all  ..near. 

.  rans}0>r*a*  io;.  vd*-n 
eff;  .‘lent  Ov  ;  ;*  lot.  ui  or.t  s, 
Network  codes  *.<  find  r.ax.ra 
flow  *  h.rou^h.  a  c.y  ,;*er  . 

"Simula*  lor:  Model  it.*-  ”  Is  inter.:*.'. 
*o  <‘i.  vmpass  arete!  -omjv  ..et.'s, 
general  simulation  j  r  oars,  ar.-i 
Simula*  1.  *.  la:  .a.  c  **xar.} 


GPSS  or  SIMSCRIPT  would  ce 
categorized  "Simulation  Modeling." 
This  category  covers  the  tools  and 
technology  of  simulation. 

H5.  Simulation  Models  "Simulation  Models"  is  intended 

to  identify  completed  models  of 
specific  systems,  however  highly 
parameterized;  e.r. ,  a  refinery 
model  or  computer  job  shop  model 
would  be  categorized  "Simulation 
Models."  This  category  covers 
the  completed  products  of  model - 
making  technology  as,  operating 
models  of  particular  object  systems. 


H6.  Critical  Path  Programs 


H7.  (Reserved.  Please  do  not  use.) 

H8.  Auxiliary  Programs  Special  purpose  utility  programs 

or  subprograms  designed  especially 
to  service  programs  in  the  above 
categories. 


H9.  Composite 


Programs  performing  combinations 
of  the  above  functions. 


INPUT  (See  also  A4)  Input  shall  contain  those  routines 

which  transfer  information  from  an 
external  medium  to  an  internal 
storage  cr  operation  area,  e.g., 
from  cards  to  an  on-line  magnetic 
tape.  Input  shall  also  include 
the  programs  which  can  perform 
either  input  or  input  and  output 
operations.  The  categories  shall 
denote  the  form  of  the  information 
to  be  inpat.  In  classifying  a 
program,  the  operation  of  input 
shall  take  precedence  over  the 
operation  of  conversion  or  internal 
information  transfer  (i.e.,  I 
rather  than  M2  or  K4). 


10.  General 

11.  Binary  Pertains  to  program  input  or 

data  input  in  the  binary  mode. 
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12.  Octal  or  Hexadecimal 

13.  Decimal 

14.  Character  (Hollerith ) 

15.  (Reserved.  Please  do 
not  use. ) 

19.  Comp,  'ite 

J.  OUTPUT  (See  also  Ak  and  K) 


JO.  General 
J1 .  Binary 

J2.  Octal  or  Hexadecimal 

J3.  Decimal 

7>i.  Chrracter  (Hollerith) 


Pertains  to  program  input  or 
data  input  in  octal  or  nexadeeimai 
mode . 

Pertains  to  program  input  or 
numeric  data  input  in  the 
decimal  mode. 

Pertains  to  program  input  or 
data  input  in  the  Character  or 
Hollerith  mode. 


A  combination  of  any  of  the 
above,  which  is  not  primarily 
one  of  the  above,  such  as  a 
general  purpose  input  program. 

Output  shall  contain  these  routines 
which  transfer  information  from 
internal  storage  to  an  external 
medium.  In  classifying  a  program, 
output  shall  take  precedence  in 
general  over  conversion.  Hie 
final  form  of  the  data  relative  to 
the  program  shall  determine  the 
category  in  which  the  program  is 
classified. 


Pertains  to  program  output  or  data 
output  in  the  binary  mode. 

Pertains  to  program  output  or  data 
output  in  the  octal  or  hexadecimal 
mode . 

Pertains  to  program  output  or  data 
output  in  the  decimal  mode. 

Pertains  to  program  output  or  data 
output  in  the  character  mode. 
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■V).  FI  ex  ling  Refers  to  routines  for  producing 

plotted  output,  either  .ia  priner 
or  via  CR'i,  or  other  special 
plotting  devices.  Routines  for 
using  plcutinfc;  devices  to  sinn.iaxe 
printing  are  also  included.  Also 
includes  routines  for  producing 
input  fcr  off-line  plotters. 

J6,  (Reserved.  Please  do 
not  use.) 

J7.  Anal op  (See  also  T6)  Refers  to  routines  which  output 

information  to  a  digital -to-analog 
converter,  other  than  that  associ¬ 
ated  directly  with  (on-  or  off-line) 
plotting  equipment,  which  will 
carry  a  -JS  classification. 

J9.  Composite  A  combination  01  any  of  the  above, 

which  is  not  primarily  one  of  the 
above,  such  as  a  general  purpose 
output  program. 

K.  INFORMATION  TRANSFER  The  category  shall  contain  routines 

whose  primary  purpose  is  the  trans¬ 
fer  of  information  from  one  internal 
medium  to  another  internal  medium 
or  from  one  external  medium  tc 
another  external  medium.  If  the 
transfer  is  external  to  internal 
or  vice  versa  the  categories  I  or 
J  should  ne  used. 

KG.  General 

Kl.  External -to-Externai  Pertains  to  the  transfer  of 

information  from  any  external 
device  to  any  other  external 
device.  This  would  be  tape -to - 
tape,  drum-to-drum,  disk-to-tape, 
etc. 

K2.  Interna  -to-Internal  Fertains  to  the  transfer  of  infor¬ 

mation  internally.  This  is  the 
same  as  relocation  of  information. 

K3.  (Reserved.  PLease  do 
not  use . ) 
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K4.  (Reserved.  Please  do 
not  use . ) 

K5«  Direct  Data  Devices  Computer-! o-compu ter  informat ion 

t rans  fer,  other  t han  i a  * :  ;e  ah o ve 
cate;  ories . 

K9.  Composite 
L.  EXECU'J  IVE  ROUTINES 
LO.  General 


LI .  Assembly 


L2.  Compiling 


L3.  Monitoring 


L4.  Preprocessing 


L5.  Disassembly  and 
De-Relativizing 

L6.  Relati vizing 


A  program  system  which  assembles 
a  program  for  a  particular  problem 
by  assigning  fixed  addresses, 
translating  from  mnemonic  codes 
to  produce  machine  language  state¬ 
ments,  possibly  incorporatin' 
library  subroutines. 

Program  systems  which  produce  a 
specific  progr  m  for  a  particular 
problem  by  translating  from  a 
problem-  or  procedure-oriented 
source  language  to  produce  a 
machine  language  program.  Inter¬ 
pretive  compilers  should  be 
included. 

Executive  or  monitor  programs, 
including  monitor  systems*  to 
control  job  sequencing,  stacking, 
log  anti  accounting  information, 
peripheral  equipment,  etc.,  with 
a  minimum  of  operator  intervention. 

Refers  specifically  to  preprocessors 
for  programs  within  the  "L"  cate¬ 
gory.  A  PERT'  (e.g.)  preprocessor 
would  be  cataloged  H8  rather  than 
L4. 
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L7.  Computer  Language  to 
Computer  Language 
Translators 


M.  DATA  HANDLING 


MO.  General 

Ml.  Conversion  and/or 
Scaling 

M2.  Sorting 
M3.  Merging 

M4.  Character  Manipulation 

M5.  Searching,  Seeking, 
Locating 

Mfi,  Heport  Generators 
( See  al  so  T1;  ) 

M9.  Composite 


This  refers  to  translation  from 
one  artificial  language  designed 
for  computing  and  data  processing 
purposes  to  anottier  such  language 
e.g.,  FORTRAN  to  PL/l.  Not  to  V.e 
used  for  translation  of  natural 
languages  such  as  English  or 
Russian. 

Data  Handling  shall  refer  to 
routines  whose  primary  purpose 
basically  is  a  mechanical  type 
of  handling  of  information  such 
as  sorting,  conversion  collating 
and  merging,  system  tape  editing 
(updating,  deleting,  adding  to,  e 


Pertains  to  any  conversion  and 
scaling  routine  (packed  or 
unpacked,  single  or  multiple 
precision)  such  as  card  image 
to  BCD,  BCD  to  card  image,  binary 
to  PCD,  BCD  to  binary,  fixed  to 
floating,  etc.  'Tie  primary 
function  of  programs  in  this 
category  must  be  conversion  ur 
scaling,  not.  input-output. 

Combined  sort. /mo rye  routines 
shall  be  included  here. 


To  be  used  for  utility  search 
subroutines.  Not  to  be  used  for 
applications  of  retrieving 
information  records  by  examining 
contents,  which  is  the  province 
of  Code  S,  Information  Retrie  ai. 
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N.  DEBUGGING 


Debugging  routines  are  those  whose 
primary  purpose  is  to  provide 
information  to  a  programmer  to 
enable  him  to  detect  errors  in 
software  or  programs.  This  will 
include  dumps  whose  purpose  is 
for  debugging,  as  opposed  to 
output,  and  search  for  debugging, 
as  opposed  to  table  look-up. 


NO.  General 

Nl.  Tracing:  Trapping, 

Snapshotting 

N2,  Dumping  (Post  mortem)  Core,  tape,  drum,  disk,  console 

printouts  (on-  or  off-line). 

N3.  Memory  Verification 
and  Searching 

N4.  Break  Point  Printing 


0.  SIMULATION/EMULATION  OF  COMPUTERS  AND  DATA  PROCESSING: 
INTERPRETERS 


00. 

General 

01. 

Off-line  Equipment 

Any  program  which  simulates 
off-line  equipment. 

02. 

(Reserved.  Please  do 
not  use . ) 

03. 

Computers 

Pertains  to  programs  which 
simulate  or  emulate  one  computer 
on  another. 

04. 

Pseudo-computers 

Simulation  of  theoretical  cr 
pseudo -computers . 

05. 

Software  simulation  of 
one  peripheral  device  on 
another 

Includes  such  programs  as  simulating 
tape  on  disk,  simulating  card  reader 
on  tape,  etc. 

09.  Composite 
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P.  EIAGNOSTICS 


Diagnostics  programs  shall  be 
ones  which  check  for  malfunction¬ 
ing  of  the  computer  or  its 
components . 


PO.  Genera] 

Q.  SERVICE  OR  HOUSEKEEPING:  These  are  routines  of  a  utilitarian 


PROGRAMMING  AIDS 

nature  which  perform  a  service  for 
the  programmer  or  operator.  Tnis 
shall  include  routines  to  find 
check  sums,  clear  memory  and/or 
compare  tapes,  etc.,  and  mechanical 
operations  involving  magnetic  capes 

QO. 

General 

Ql. 

Clear/Reset  Programs 

Q2. 

Check  Sum  Accumulation 
and  Correction 

Q3. 

Rewind,  Tape  Mark,  etc. 
Programs 

Q4. 

internal  Housekeeping; 
Sace,  Restore,  etc. 

05. 

(Reserved.  Please  do 
no  1.  use .  ) 

06. 

Program  Documentation: 
Flow  Charting,  Document 
Standardization,  etc. 

NON- 

-NUMERICAL,  MATHEMATICS 

Logical  functions,  logical  oper- 

ations,  logical  calculi  and 
algebras,  symbol  manipulation 
and  manipulation  of  non-numeric 
quantities. 


R0. 

General 

HI. 

Formal  Logic 

Rf>. 

Symbol  Manipulation 

Includes  comiinational  generators, 

permutations,  etc. 


R3.  1. 1st  and  String  Processing 


S.  INFORMATION  RETRIEVAL 


Programs,  or  systems,  for  cataloging 
data,  all  related  to  one  field  of 
interest,  so  it  can  be  displayed, 
with  accuracy  and  speed,  on  an 
output  unit  at  any  time  it  is 
needed.  May  or  may  not  be  real¬ 
time.  Examples  are  classification, 
indexing  and  machine  searching 
systems. 


SO.  General 
T.  APPLICATIONS  AND 


Application-oriented  programs 

TO. 

General 

TL. 

Physics  (including 
nuclear) 

T2. 

Chemistry 

T3. 

Other  Physical  Sciences 

T4. 

Engineering 

T5. 

Business  Datei  Processing 
(See  also  H  end  M6) 

T6. 

Manufacturing  (non-data) 
Processing  and  Process 
control  (See  also  J7) 

'17. 

Mathematics  and  Applied 
Mathematics 

T8. 

Social  and  Behavioral 
Sciences  and  Psychology 

T9. 

Biological  Sciences 

Programs  or  subroutines  which 
solve  specific  problems  in  fields 
other  than  computers. 


Programs  pertaining  to  physical 
sciences  not  listed  above. 
(Geology,  Astronomy,  etc.) 
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U.  LINGUISTICS  AND  LANGUAGES  Programs  which  work  with  natural 

languages,  usually  translating 
from  one  language  to  another 
language.  The  rules  of  these 
languages  reflect  and  describe 
current  usage  rather  than  prescri 
usage.  Example  would  be  translat. 
English  to  French. 

UO.  General 

V.  GENERAL  PURPOSE  UTILITY  SUBROUTINES 

VG.  General 

VI.  (Reserved.  Please  do 
not  use. ) 

V2.  (Reserved.  Please  do 
not  use.  ) 

X.  DATA  REDUCTION  Many  laboratory  or  field  tests 

and  experiments  automat  leal Ly 
record  data  either  at  the  site 
or  by  way  of  telemetry.  Programs 
in  this  category  will  accept  sun 
digital  data  and  perform  the 
necessary  functions  of  deoommita- 
tion,  scaling,  ca Li  orating, 
evaluating  and  test  analysis. 

Some  of  the  programs,  especially 
X4,  Xcj  and  X1',  might  t-e  predom¬ 
inantly  of  type  P,  E,  F,  G,  and 
could  possibly  i>e  found  under* 
those  categories,  "he  rtvvrait.r 
are  either  for  post -pro  *ecs in 
or  they  may  be  on-line  <-ra*  1 .  • 
in  real  time. 

XO.  G.-tvru! 

X:.  K*i\  rmat ♦ 5 no.  Procram  separates  the  arias, ies 

Decommu  t  at  ion,  a:  id  often  \  er*s  them  o  Vvrr.» 

Error  ria(*nosis  puter  words  or  higher  le  el 

language  variables.  Miss  ire  or 
erroneo us  data  are  identifies. 
Output  consists  of  ordered  data 
and  otitinc  i r. format ic:.. 


Ij-lh 


X2.  Editing 


X3.  Calibration 


XL.  Evaluation 


X5.  Analysis 


Xu.  Simulation 


US EH  ASSIGNED 


ALL  t>TKERS 


Output  from  category  XI  is  used. 
Bad  data  Eire  cast  cut,  tr.issi re¬ 
values  inserted,  wrone  alues 
corrected.  Process  is  ei’her 
automatic  or  by  parameter  aris, 
or  both.  Output  is  called  .lean 
raw  data. 

Data  are  scaled  linearly,  Then 
calibrated  to  obtain  fun -tier, 
values  in  physical  units. 

Output  is  called  "clean  cal  it rateu 

data. " 

All  necessary  computation  on  the 
data  is  performed  to  present  them 
in  a  form  suitable  for  engineer! i.g 
or  scientific  evaluation. 

All  computations  necessary  to 
analyze  the  outcome  of  the  tes* 
or  experiment.  Also  referred  t> 
as  Time  Series  analysis. 

Programs  which  generate  artificial 
data  to  be  used  as  the  thee re ti  *ai 
test  model  or  to  te  used  for 
checkout  of  programs. 

Ill  is  category  is  for  the  use  of 
individual  installations  who 
have  groups  of  routine-,  which 
do  not  fit  well  under  ary  of 
the  above  categories.  5  -urines 
from  this  -category  must  ■ 
re-classified  on  an  indi.iuual 
t  asis  when  and  if  tie;,  ur  ■  u.r- 
mitte  i  to  the  User  In  !  i  •  rarv , 
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ADDENDUM  A 


1.0  IDENTIFICATION 

1.1  Title 

Real  Zeros  of  a  Single-Valued  Function 

1.2  Identification  Name 
DO-NKL - FCNZERO 

1.3  Classification  Code 

DO  -  Operations  on  Functions  and  Solutions  of  Differential 
Equations,  General 

1.4  RCC  Identification  Number 
D0001R00 

1.5  Bntry  Points 
FCNZER0 

1 .6  Programming  Language 

Language :  FORTRAN 
Routine  Type:  Subroutine 
Operating  System:  DRUM  SCOPE  2.0 

1.7  Computer  and  Configuration 
CDC  3800 

1 .8  Contributor  or  Programmer 

Janet  P.  Mason,  Code  781 3 >  Research  Computation  Center, 

Mathematics  and  Information  Sciences  Division 

Harold  L.  Toothman,  Code  4319H,  Naval  Analysis  Staff, 

Office  of  Director  of  Research 

1.9  Contributing  Organization 

NRL  -  Naval  Research  Laboratory  -  Washington,  D.C.  20390 


Prucra::,  A.  ali  as  i  I  i  ty 


1.10.1  Su-mit  tal ;  Program  write-up,  Fortran  source  deck,  source 

listing 

i  .1C. 2  Os.  File:  RCC  Prr  -ram  Library 

1.11  Verifi cat ion 

The  fol sowing  two  problems  were  used  as  test  cases: 

( i )  Seventy  zeros  of  the  function  sin(l/X).  See  Section  7-0. 
(£)  F(e/v?)  _  (Rq*  -  Rr2)T!  =  Q 

Ho/'h)  (%  -  FTi)Tt 

where  Uiei6l60  and  Tg,  RQ1 ,  RQ2,  Rrl,  and  Rf2 
are  known. 

1 . 1 2  Date 

1  Aug"  is-  1968 

ruw  oss 

d .  1  Do script) on  of  the  Routine 

To  find  the  real  zeros  of  a  single-valued  function  of  one  real 
variable  sy  a  modified  Method  of  False  Position. 

Find  wi thiii  a  specified  tolerance  &  those  points  along  a  given 
closed  interval  [u,b]  for  which  the  functional  value  f(x^) 
satisfies :  |f(*.  il<«.  isach  point  and  the  corresponding 

in  a  table: 

f  (  )  >  •  •  •  >  X]QJ>  1  (Xn) 

wiiere:  a;£xr,<xn£b 

P  .2  Fronl em  background 
See  the  authors. 


functional,  value  are  stored 
X0,  t(x0),  XL, 


USAGE 

3 . 1  Calling  Sequence  or  Operational  Procedure 

FCNZ11R0  ( F0X,  X0WER,  UPPER,  T0LE,  MAXRTS,  NUMB,  ANS) 


A-? 


F0X  is  the  name  of  a  FORTRAN  function,  to  be  written  by  *.r.e  ..ser, 
which  evaluates  the  function  whose  zeros  are  ic  be  calculates. 
F0X  must  be  declared  external. 

X0WER  is  the  lower  bound  of  ne  given  interval. 

UPPER  is  the  ’typer  bound  of  the  given  interval. 

i0LE  is  the  tolerance  set  upon  both  x.  and  f(x..  ). 

MAXRTS  is  the  number  of  zeros  (less  than  513)  the  user  expects  '  e 
find  in  the  given  interval.  The  user  van  make  MIX?. Tv  -:q;a: 
to  0,  thereby  leaving  the  determination  of  the  number  of  zeros 
( e .  g . ,  roots)  entirely  up  to  the  routine. 

NUMB  will  contain  the  number  of  zeros  actually  found.  I.'JMr  r.av  or 
may  not  equal  MAXRTS. 

ANS  is  the  name  cr  +he  table  of  answers.  The  user  should 
dimension  ANS(2050)  or  ANS(h  x  MAXRTS). 

j . 3  Space  Required  (Decimal  and  Octal ) 

3.3.1  Unique  Storage:  1020  octal  (528  decimal)  locations 

exclusive  of  computer  system  library  Surety 


j.3.2  Common  Plo  eks : 

C0MM0N/A,/  INPUT,  0UTPUT,LIM, VAL , Z0 , Zl ,  FZO ,  FZl 
INPUT  \ 

OUTPUT)  LIM  -  integer 

VAL  l 

Z0  V-  real 

Zl 

FZO  I 

FZl  J 

length:  10  •  ctal  (8  dec)  locations 
3.3.3  Temporary  Storage:  none 
3.*+  Messages  and  Instructions  to  the  Operator 


Error  Returns.  Messages,  and  Codes 


.(>  Informal i v e  Messages  to  the  User 


firdinar L  1  y  there  is  no  printout.  If,  ..ovever,  ihe  limi.  (l  o'1  on 
the  number  of  iterations  in  the  False  Position  method  has  eor 


reached,  or  if  the  two  points  xg  and  bounding  the  trie  zero 

sat isfy: 

J(xs  "  0.0000000001, 

then  the  following  information  is  printed: 

nut  ter  of  iterations,  Xg,  f(xs),  Xt,  f(xt). 

The  last  point  calculated  is  taken  to  be  the  zero  (e.g.,  root) 
and  is  stored,  as  is  its  functional  value. 

3.7  Inpu* 

N..  data  input. 

1.8  0  n  nut 

— .  - 

Ordinarily,  none;  see  3*c. 
j.9  Format 

not  applicable 


3.10  External  Routines  and  Symbols 


FCNZHR/  -  Q8QDICT. 

CONVERGE 


(system  library) 
(deck ) 


CONVERGE 


THEND. 

Q8QDICT. 

STH. 

QNSINGL.  t 


(system  library) 


3.11  Timing 


1.8  seconds  for  sample  program  (e.g.  to  find  at  least  70  zeros  of 
the  function  sin(l/j£  .00001  5  x  i  1,0);  see  5.0 


3.12  Accuracy  -  -  - 

3.13  Cautions  to  Users 


If  the  programmer  wishes  to  find  the  first  n  zeros  in  a  given 
internal  he  should  let  MAXRTS  =  0  (see  3*2)  in  his  subroutine 
call  and  pick  out  the  first  n  zeros  from  the  table  of  answers. 

If  MAXR'IS  were  to  be  made  equal  to  n  (and  n  is  less  than  the 
actual  number  of  zeros  in  the  given  interval),  ii  is  quite  likely 
that  the  zeros  found  would  not  be  the  first  n  zeros  in  the  interval. 
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3.1 k  Program  Deck  Structure 
J0i:  card 

FTN  card  *  __ 

main  program  deck  (contains  call  to  FCNZER0)  {( ^ee 

Function  deck  (a  single- valued  function  of  one  real  '.aria-  i  ’to: 

subroutine  FCNZER0  l  DC-NKL-FCNZER0  (See  Section  5.l) 
subroutine  CONVERGE J 

SC0PE  card 
L0AD  card 
RUN  card 

3.15  References  -  Literature  -  Appendices 

Introduction  to  Numerical  Analysis  by  F.  E,  Hildebrand 
METHOD  OR  ALGORITHM 

After  the  number  and  approximate  location  of  zeros  are  determined  ij 
successive  sectioning,  the  function  employs  a  modified  Method  cl  False 
Position  to  determine  the  exact  location,  within  the  giver,  tolerance, 
of  these  zeros.  False  Position  (iterated  linear  inverse  interpolation 
ai though  known  to  converge,  may  do  so  very  slowly.  To  hasten  con  r- 
gt-hoo,  the  usual  False  Position  method  has  been  modified;  the  stops.— 
in  interpolation  is  adjusted  whenever  it  becomes  relatively  smal .  and 
t.'ie  last  two  successive  approximations  of  the  zero  have  the  same  stgn. 
Fonvergence  is  satisfied  when  either 


(0  |f(x.  )\<S  ft  Jxt  -  ^ 


where  f(x  )<  0  and  f(x  )  ,>  0,  or 

S  T» 


(£)  f(x.)  «  0. 


no  no  ono  n  on 


5.0  FLOWCHART  AND/ OR  SOURCE  LAKQUAOE  LISTTNO 


C 

C 

C 

c 

c 


SUBROUTINE  FCN28P0(LC.L0.U3.TL#NR,SC,LS  ) 

'DENT  nUmBER  -  DCOOlROO 

TITLE  -  REAL  ZER6S  0F  A  S I vGlE. VALUED  F  NCTIQN 
IDENT  NAME  -  D04NRL-FCN7ER3 
LANGUAGE  *  FORTRAN 
CGMPUTER  -  CDC-36  00 

CENTR I BLT6R  -  „ANET  P.  MAS3N,  C0DE  7813.  RESEARCH  C0MPUTAT 10N 

CENTER,  MiS  DIVIS10N 

ih  AF  OLD  L.  T03THMAN,  C0D6  4319H,  NAVAL  ANALYSIS 

STAFF.  0FFICE  0F  DIRECT8R  0F 
RESEA-  CH 

ORGANIZATION  -  NFL  -  NAVAL  RESEARCH  LAB  RAT6RY  -  |>  ASH  1 NG  T0N  •  D.C. 

2  039o 

DATE  -  1  AUGUST  1968 

PURPOSE  -  T0  FINE  THE  REAL  ZER0S  0F  A  S : NGlE- V ALyE D  FUNCTION  0F 

ONE  real  VARIaELE  by  A  ‘8DIFIFD  METHOD  F  FALSE  F»€  S I T 1 0N 

DIMENSION  BND ( 4  ) ,  LS T ( 2 050 1 

type  real  lbout,  input, lb. lst.m.lc 

type  integer  tf .fsc.r.s.t ,V,SC 

CGMP0N/A/ INPUT, 0L TPUT, L IM. vAL.Zo.Zl.FZo -FZi 

LIM=iooO 

DELTaXiO.O 

V  A  L  =  0 . 02 

K4s0 . 367879441 

Tfc=FSCsSC*R=V=C 

DxsCEL  TAX 

IF(DX. EC. 0)1, 2 

1  DXSUB-LB 

2  1F(NR.EC.0)3,4 

3  NP=5i3 

4  DLMMYsINPCTslB 
LB0L  T  =  LC ( DUMMY  ) 

5  SC=R=0 
8ND<1)=L B 
'■iND  (  2  )  *L  B€UT 

6  IF(BNC<1 ) .EO.UE 17 .8 

7  I F ( A0S ( PNC ( 2 ) ) . L  T .TLJ9.10 
9  LSTIR+i )«BND<1 > 

LST (R*2)eBND(?) 

SCsRC+1 

10  GG  TP  20 

8  I N  PL  T  sBND ( 1 ) »DX 

!F( INPUT. GT.UBJH  .12 

11  IK  PUT  sup 

12  0ND(3).ruPHY3lKRLT 

Bn  D ( 4 ) sl  C ( Dummy > 

IF ( BNC ( ?  )  .EG. 0)17 .13 

13  IF ( 0  N  T ( 2  >  •LT,0>}4 ,15 

14  IF  (  P  NJ  T  (4  ) .  G  T  ,  0  )  1'  7  *16 

15  IF(0Nr(4).LT.())l'7,16 

16  GG  TO  i 9 

17  DG  10  1*1.4 
LST (R  +  i ) *END ( I  ) 

18  cgntinuf 


10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 
30 
40 
50 
60 
70 
80 
90 

100 

110 

120 

130 

140 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 
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RsR*4 

SC*SC+1 

IF1SC. Gf .513)25.19 
1’  BN0(i)sPNC(3) 

BN  D  f  2  )  =BNC ( 4 ) 

GG  T 0  6 

20  If(SC.EQ.FSC)?l ie2 

21  Th*TW41 

^  GG  T P  23 

fsc=sc 

GG  Tp  5 

25  Ts ( SC-1 ) *4*1 
oe  32  S«1.T.4 

IF<LST(S41)-E0.o)26.27 

26  LST( v*i )3LST(S) 

LST  ( V*2)3LST(S'*]  ) 

GG  TP  31 

27  IF  C I  ST(S  +  1)  .'_T. 0)28, 29 

28  Zo=!Sns)J  G  ?[  =LS  T  (  s*i  ) 

%m\VlY'i9  G?1  =LST  (  S  +  3 ) 

2’  21*1  ST(S)S  F7.l«LST(S4l) 
70*LST(S42)S  G7C=IST(S*3) 

30  CALL  C8N VERGE ( L C  #  TL  ) 

L  S  T  (  V  ♦  j ) *  I NPU  T 

LST ( V+2 ) sCUTPuT 

31  ViV*2 

32  CGNTINuF 

33  end 


400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

340 

530 

560 

370 

380 

590 

600 

*10 

620 

630 

640 

650 

660 

670 

680 

690 

700 

710 
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s t  b v p i  t r n e  o:\v?*ge(lcn,tl v>  eoo 

TYPF  RE*l  K3  ,  I  Ml  T.lCTv  810 

CFMvpi\/A/lNP'  T,  'L iPlT.LJK,  AL./0.71.FZU.F71  820 

FACTPRsl.O  830 

K2  =  0  840 

Ff  PR|;T*k3«1. 0^1  Ft  850 

34  0LMk'Y=^PLT=(F2'.»Z0-F/0*Zl>/(F71-FZO)  860 

0 L  T P !j T s l  C M  P P F  v  )  870 

K?=K?+1  880 

1 F  <  *pc ( PUTPU  t ) . *  l.TLN.AfD.ARS(71»70).LT.TLN)50,35  890 

35  RAT IPrAPSf { F  « P I N -  INPUT )/ (Zl-Zp  )  )  9oO 

!F(RATIP,LT,VAL.AND.I:f:Pi-UT*0UTPUT.GT,o)36.42  910 

36  FACTPRsFACT0R»1  r  .  f|  920 

37  IP{RATIP*cFACT€-ii  .q)  .GP,  1.0)38, 39  930 

38  FACTPP:FACT0P/2  .  C  940 

GF  TF  37  950 

39  SFECMMlT  960 

AlSPsSL-TPIT  970 

DLMYY-fKPLT=[NP'  1 -F AC  I 0t* C  r  PR  r 1  -INPUT )  980 

0t.TPtir=ic^  <p.jm*y )  990 

K?  =  K2+1  1000 

IF  (ALSF.LT . n ) 4  C  <  a l  1010 

4 0  FfRINsZf'sSPFC  1020 

FFRF|'TrFZPsALeG  1030 

GP  TP  A?  1040 

41  FCR’a sZI =SPFC  1050 

FFRFi  T-FZ1*ALRG  1060 

42  IF (  F  U  7  P(  T.EQ.0)5C.43  1070 

A3  IFfPCTPI  T.LT.nM*  ,45  1080 

44  FFRl\sZr*INPi.T  1090 

FFRFijTsFZOap'.TF.  T  1100 

GF  TP  4A  1110 

45  FpRI\r71sINP_T  1120 

FERFI;T!FZl*0tTF,  1  1130 

46  IF  (  Ape  (PUTPUT  )  .1  1  .  TL  hi.  AA'D.  APS(Zl-ZO)  .LT.Tl  N)50.47  1140 

47  IF (K2 . GT . L IM. CP ■ AfiS (( 7q-Zi )/Zn ) .LT .0 .0000000001 >48 .34  1150 

48  PRIM  4?,K2,Zn.f20*Zl.F?l  1160 

49  FGRYATdX, I4.3X .E 1 7 . 10.3X. :17.10*3X.F17.10#3Y.E17.10>  1170 

50  FND  ll8o 


6.0  COMPARISON 


No  comparisons  have  been  made  with  any  other  programs. 
7.0  TEST  METHOD  AND  RESULTS 


The  main  program  (’TRIAL)  and  function  evaluator  (FOX)  used  to 
evaluate  the  seventy  zeros  of  the  function  sin(l/x),  along  with 
the  results  as  printed  out  by  the  main  program,  follow  this  page. 


MAIN  PROGRAM 
(CALLS  FCN2ER8 J 


program  trial 

001 

External  fox 

S001 

DIMENSION  ANS< 2050 ) 

T001 

NT*0 

U001 

Fermat  <inu 

V001 

PI»3. 1415926536 

002 

maxrts«70 

003 

T6lE*5 , 0G»4 

004 

XpwPRel . 06-5 

00? 

UPPER*1.0 

006 

CALL  FCN2ERB<Fex,X0WER. UPPER, T0lB,MAXRTS.NUMB.ANS) 

007 

NLM0=NUMB*NUM9 

006 

PRINT  3 

009 

format ( ini, i7x, *x». 15x,*sin( i/x)*, Ox, •( i/x  /pi •//> 

S009 

De  1  j«i , NUMB, 2 

Oil 

FINAN*ANS< J) 

012 

REC*1.0/(FINAN*PI ) 

013 

fun*ans< J*l> 

014 

PRINT  2.  F  INAN.FLN.REC 

015 

F0RMAT(OX,El7.1o,3X.El7.lO,3X.F9.3) 

016 

NT*nT*i 

017 

IF(NT.0E. 56)5,1 

018 

PRINT  4 

019 

NT*0 

010 

continue 

021 

END 

022 

function  Evaluator 

(EVALUATES  SIN  1/X) 


function  rox(x) 

01 

x*i, 0/x 

02 

F€X«SIN(X) 

03 

END 

04 

A- 10 


X 


S  I  M  ( 1  /  X  ) 


<1  /*)/(■  I 


1  ..1340732483-004 
*3  .  012754^73-004 

8  .0365487938-004 
7 .  C8S3C61138-C04 

9  .41745  2^791  -004 
1  .  1090931983-003 
1  .1449995o78-C03 
1  .?48274r6?i-0Q3 
1  .4336280048-003 

1  .97708 0O34e-OO3 

2  .0669473475-003 
2.1080124913-003 
2 . 2736420  442-003 
2.357851OO57-C03 
2.8420519696-003 

3  . 09O3e99<.90-003 
3 .2460574004-003 
3 .49791  04j55-003 
3 .417158^596-003 
3 .835058X  762*003 
4.0808561100-003 
4 .8970751  Hi 0-C03 
5.l340?8xO3e-C03 
5 .39t'062xoo3-C03 
5.4880962515-003 
5 . 7874524493-003 
5 . 89462  7x903-003 

6  .  no5847r  O26-C03 
6 . 121  343Q39Q-0P3 
6 . 2413703614-003 
6.36619/7443-003 
6.49611759x7-003 
6.6314559734-003 
6 . 7/2543931  4-  003 
4.91977V7/4J-C03 

7  .  07?559*{)i  C-C03 
7 .  7343156302-003 
7  .  4 c  7571 7666- C03 
’ . *  7880  920X6- C0O 
7. 7636557472-003 
7.95774  76522- C03 
8. 161 7934 0 15*003 
* .376575X450-003 
8.6O29699H4-C03 
8.8419412966-003 
9. 89456 81944-COJ 
9.3e70*54«i2-C03 
9.6457531932-003 
9.C4714J5430-0O3 
1 . 0268057192-002 
1 . 061 0335341-002 
1 . 0976211530-00? 

1 . 13682X0386-002 
1 .1789239332  002 
1 . 7P426O7491-C02 
1 . 27323VXSX2-C02 


2.1490885245-004 
-1,0243733414-005 
1  .5Q77646821-005 
-1 .8468816311-004 
-8  .1  76276*  255-005 
-3.9168868260-004 
2 .1455052  774-006 
-8  .3513441496-007 
1 .1254008835-004 
-9 . 0987305157 -007 
-6 .9291007260-006 
1  .6475331383-008 
-9.5486197339-009 
1.0244172534-005 
7.0021836961-005 
7 .0988605786-004 
2.4481041811-004 
-1 .6478938050-005 
-2.0050393731-005 
4.4779409874-007 
-8.0076595789-006 
3.7483079580-Q07 
6.9748522890-005 
1.4991408/01-006 
1 .0738559629-004 
-7.9983146861-007 
-1 . 9697789125-006 
4.2551280952-006 
6.9400994107-007 
1.6471629518-006 
-5, 0897(335038-007 
-5.9965619547-005 
-2.5541777943-007 
-1 .4908969697-004 
7.5194693636-006 
1 . 3659880128-004 
-6.7215552774-007 
2,9698348773-004 
-4,1691094161-005 
- 8 . 5165956989-006 
-7.8577457482-006 
2. 1738060983-005 
1.5711702872-006 
1.9633700360-007 
-1 . 7573149164-007 
2. 365932*960*007 
-5.9985557916-008 
-1  .  O03l9i>,735-005 
4 . 0830466925-004 
-3. 464846*851-005 
-6,1 >37443863-005 
7 . 1035496832-005 
-1 . 2774274897-006 

-1 . 1304809571-004 

2 . 9267685022-006 
x . 31 9Q64631 2-CC  7 


2386 . 0  C  0 
635  .  OC  0 
632 .  OCO 
449. OCO 
338  *0f 0 
2 87. OCO 
27b. OCO 
255. OCO 
222. OCO 
161. OCO 
154. OCO 
151. OCO 
1<U.  OCO 
135. OCO 
112. OCO 
103.  OCO 
98- OCO 
91.  OCO 
88. OCO 
83. OCO 
78. OCO 
65. OCO 
62. OCO 
59. OCO 
58. OCO 
55.  OCO 
54. OCO 
53. OCO 
52. OCO 
51.  OCO 
50. OCO 
49. OCO 
4b. OCO 
47. OCO 
46. OCO 
45. OCO 
44. QtO 
43. OCO 
42. OCO 
41. OCO 
40. OCO 
39. OCO 
38. OCO 
37. OCO 
36. OCO 
35. OCO 
34.CC0 
33. OCO 
32. OCO 
31. OCO 
30. Of  0 
29. OCO 
26. OCO 
27. OCO 
26. OCO 
25. OCO 
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1  .  3 ?<?•  2 ^  (J4V1  3-C02 

3. 9857666696*005 

24 . 0C  0 

1  ,  3e3  95s6n4?g  ->C02 

8. 2978187350-007 

23.0C0 

1 . 44*8*31 136-002 

2.5995541364-007 

22. 0C0 

1  .5i57ti5567-C02 

8.4384544773-006 

21. 0CQ 

1 , 591  5494  7?5-002 

-1.6415459126-006 

20.0C0 

1 ,6753(44592-002 

-3.8234935142-004 

19.0C0 

* 

1 . 7*P3fo?7i 1-C02 

-4.6319002506-007 

i8.oro 

1 ,e7?4127098-C02 

4.6053147162-005 

17.0C0 

1 

1  .9«<J436«aie-C02 

-3.3347425272-007 

16,0(0 

i 

2  , 1??Qt0<3391-C02 

-1.1034033355-004 

15. 0C0 

■? 

2 , 27366 14639-002 

-3.7546624037-004 

14.0C0 

j 

2.4485578572-002 

3.3811357182-004 

13. 0C0 

1 

2.6525637247-002 

-1.9042347792-005 

12.0C0 

i 

?.P93740Oo«5-002 

1.7591417583-004 

11. oco 

4 

3.1831C2205S-C02 

-3.3004473891-005 

10.0C0 

1 

3.5367671674-002 

8.5173788376-005 

9.0(0 

3 . 978863no7l-C02 

-5.9563698214-005 

8.  OCO 

j 

4.5472924465-002 

4.0432819156-005 

oco 

5.3051686926-002 

-1,4648542506-005 

6.0(0 

6. 3662016117-002 

9.5932071684-006 

5.0(0 

7.0577506*56-002 

-5.576G636481-006 

4.0(0 

1.061 0332581-001 

2.7021185983-006 

3.0(0 

1 .591549*642-001 

-9.209718*435-007 

2.0(0 

3 . 183058Q482-C01 

8.5221142264-008 

1.0(0 

8.0  ’'BUMP 
Mom 
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ADDENDUM  B 


PROGRAM  BRIEF 

TITLE.  Real  Zeros  of  a  Single-Valued  Function 

IDENTIFICATION  NAME.  DO-NRL-FCNZERO 

PROGRAMMING  LANGUAGE.  36OO  FORTRAN 

COMPUTER  AND  CONFIGURATION.  CDC-38OO 

PROGRAM4ER.  Janet  P.  Mason,  Code  7813 

Harold  L.  Toothman,  Code  U319H 

NRL  ORGANIZATION.  Mathematics  &  Information  Sciences  Division 

Research  Computation  Center 
Code  7813 

• 

DATE.  15  January  1968 

PURPOSE,  To  find  the  real  zeros  of  a  single-valued  function  of  one  real 
variable  by  a  modified  Method  of  False  Position. 


PROGRAM  STATUS.  Complete 


LIBRARY 


ADDENDUM  C 

PROGRAM  ABSTRACT  COVER  SHEET 


(SUBMITTAL  FORK) 


catalog  identification  do-nrl-fc?ki~ro 


REVISION 


TITLE  (DESCRIPTIVE)  : 


SOURCE  LANGUAGE (S) 
AND  DIALECT(S)  : 


OPERATING  SYSTEM 
AND  VERSION  : 


DRUM  SCOPE  Ver 


COMPUTER  &  CONFIGURATION  :  CDC-3800 


CONTRIBUTING  ORGANIZATION  :  ."RL  -  Naval  Research.  labor  at  or 

Janet  P.  Mason,  Code  7813 

SUBMITTOR/ PROG  RAMMER  :  H.  L.  Toothman.  Code  REVISOR  : 


exclusive  of  computer 


DATE  WRITTEN  : 


ORIGINAL  IS  January  '68  REVISION 


PROGRAM  MATERIAL  :  WRITE-UP  PAGE  COUNT  : 


SEQUENCED  SOURCE -CARD  COUNT  : 
M/T  RECORD  COUNT  :  _ 


P/T  RECORD  COUNT  : 


OTHER  AVAILABLE  DOC!  (ENTATION  : 


ORIGINAL  /  REVISED  PROGRAM  CATALOG  ABSTRACT 


A  subroutine  to  find  he  real  zeros  of  a  single-valued  function  of  one 
real  variable  by  a  modified  Method  of  False  Position. 

Entry  point:  FCFEERO 


For  Completion  by  Distribution  Agency  : 
Original /Previous  Revision  :  UGN 
Lacest  Revision  :  UGN 


Accessed 

Accessed 


C-l 


DOCUMENT  CONTROL  DATA  R  &  D 


. . .  ,./  f,./--  I,,- 1,  ,l 


Naval  Research  Laboratory 
Washington,  D.C.  20300 


'  li  vnten-J  t*  *»»  n  tin  :  >  v  ■  •  f  «i  /  /  r**porf  /  .  <  f.ivw/infi 


1«.  '<1  CCi.  f  M  I  >  F4  :  ’  v  f  A  i  i  i  ■  A  I  it'!, 

UNCLASSIFIED 


RESEARCH  COMPUTATION  CENTER  PROGRAM  PUBLICATION  GUIDE 

OESCRIPUVE  NOTES  (Type  ol  repofl  and  inc/u.i  IVP  dates  ) 
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